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Abstract. We study the semiclassical propagation of squeezed Gaufiian states. 
We do so by considering the propagation theorem introduced by Combes- 
cure and Robert [T] approximating the evolution generated by the Weyl- 
quantization of symbols H. We examine the particular case when the Hessian 
H" (Xt) evaluated at the corresponding solution Xt of Hamilton's equations 
of motion is periodic in time. Under this assumption, we show that the width 
of the wave packet can remain small up to the Ehrenfest time. We also deter- 
mine conditions for "classical revivals" in that case. More generally, we may 
define recurrences of the initial width. Some of these results include the case 
of unbounded classical motion. In the classically unstable case we recover an 
exponential spreading of the wave packet as in pp. 



1. Introduction 

Localization of quantum states in phase space is a prerequisite in some semiclas- 
sical treatments of quantum evolution. In the classically chaotic case, the width of 
an initially localized Gaufiian increases exponentially [I] up to the so called Ehren- 
fest time Te, i.e., the time up to which quantum dynamics can be approximated 
by classical dynamics. In the case of regular classical motion it can be shown that 
this width grows algebraically in the semiclassical parameter h up to Te. 

In certain applications, it is necessary to propagate states semiclassically for 
long times. This demands control of the width of the state. An example is the 
construction of quasi-modes proposed by Paul and Uribe [2]. One refers to a quasi- 
mode as a state ip which is a solution of a corresponding spectral problem of an 
operator H up to some small discrepancy 5, i.e., 

Hip-Eil) < 5 

which insures, in the case of discrete spectrum, that there exist at least one eigen- 
value of H in the interval [E — S, E + 5}. The approximation proposed by Paul 
and Uribe uses the semiclassical propagation of coherent states over closed classical 
trajectories leading to the well known Bohr-Sommerfeld quantization rule in one 
dimension. The eigenvalues E n of the Hamiltonian H are given by the quantization 
condition 

(1) J pdq = 2<!rh(n+^J + O (h 2 ) , n e IN, 

where the energy shells % = {(p, q) : H (p, q) = E} are closed curves. 

Certain examples suggest that there exist systems (other than the harmonic 
oscillator) where the propagated width of an initial Gaufiian remains small for long 
times. Such a behavior is exhibited by the propagation of Gaufiians generated by 
some perturbed periodic Schrodinger operators like the Wannier-Stark Hamiltonian 

h 2 

#ws (e) = ~Y Ax + Vr ^ + ex 
in the limit of small perturbations e. Here Vr (x) denotes a periodic potential with 
respect to a lattice Y = ifi . It is known that the band structure of the spectrum of 
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the unperturbed operator is preserved for small enough perturbations e. Numerical 
studies [HI H] show that an initially localized Gaufiian in momentum space defined 
on such an energy band apparently remains Gaufiian for long times. The evolving 
states carry out oscillations. The center of the Gaufiian can oscillate in position 
space describing so called Bloch oscillations [§]. Alternatively, the width of the 
Gaussian can oscillate, whereupon the center remains fixed, describing so called 
breathing modes. In both cases, the state returns to the initial state after an 
oscillation period up to a small error. 

We study the evolution generated by a class of operators on L 2 (R d ) defined as 
Weyl quantizations of classical symbols H (X) with the property 

(2) H"{X t ) = H"{X t+T ), VteR, 

where X t denotes the solutions to Hamilton's equations of motion and H" (X) 
is the Hessian of H with respect to X. In one dimension, the condition ([2]) is 
satisfied by bounded classical motion or unbounded motion in a periodic potential. 
We particularly study the evolution of initial Gaufiian (or squeezed) states semi- 
classically, i.e., asymptotically as h \ 0, when t /" oo. We focus our attention to 
the spreading of such wave packets. 

In sections [l.ll and [L2l we shortly review the semiclassical propagation of Gaufiian 
coherent states. In section [T73l we give some known results on the validity of the 
approximation. We then make, in section |2~T1 and section statements about the 
approximate Gaufiian state given by this semiclassical propagation using Floquet 
theory. These properties are then brought back to the true evolution in section [31 

1.1. Preliminaries: semiclassical propagation of wave packets. We will 
work in the context of self adjoint operators defined on L 2 (R d ) that are K— Weyl 
quantizations of symbols. To a smooth (C°°) classical symbol b(X), i.e., a func- 
tion on the phase space T*R d = R 2d , there corresponds an operator on L 2 (R d ), 
b := Op^ [b], defined by 

(2irh) Jrzi V 2 / 
The following conditions will be assumed : 

(1) The classical Hamiltonian H : H 2d — > R is a smooth function. 

(2) H £ S (to), i.e., for all multi-indices a, (3, there exists K a ^ > such that 

\d^ q B (p, g) \<K a ,p(l+ \p\ 2 + \q\ 2 ) ~ . 

(3) The corresponding classical equation of motion is given by 

where J is the symplectic unity 

7 - ( ~ ldxd ^ 
J ' V ^xd J ' 

and H' is the gradient of H with respect to X. Furthermore, we denote by 
$^ : R 2d — > R 2d , X ^ X t — & H (X Q ) the corresponding classical flow. 

(4) The h— Weyl quantization of H (X), H :— Op^ [H], is an essentially self- 
adjoint operator on L 2 (R d ) and generates a unitary time evolution Vt G R, 

U (t) : L 2 (R rf ) -> L 2 (R d ) 
4 (0) » 1> (t) 
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corresponding to the Schrodinger equation 
(3) ih d -± = Hip, 



i.e., we will write U (t) = e « 

The Weyl calculus also allows a representation of quantum mechanical wave 
functions on phase space. This is given by the Wigner function of the state u G 
L 2 (R d ) , 

(4) W[u](p,q) := J u (q + |) u (q - |) dy. 

The Wigner function of a Gaufiian defines a positive measure on phase space. 

1.2. A semiclassical propagation theorem. With these assumptions we give 
a short summary of the method of semiclassical propagation of coherent states 
introduced by Combescure and Robert [I]. Similar constructions have also been 
considered in the past. See, e.g., Hagedorn [6, 7\ and references therein. 

The idea is to expand the exact Hamiltonian H along the classical flow generated 
by the symbol H up to second order, whereupon the approximate time dependent 
Hamiltonian, iJ 2 (t), is the H— Weyl quantization of 

H 2 (t,Y) := H(X t ) + {Y-X t ) T H'(X t ) + ^(Y-X t ) T H"(X t )(Y-X t ). 

The propagation of normalized wave functions (squeezed states of the form 



(5]P<i-^(x) := dct ^(f))% Mp r (^)+(^) r 4(^)),ZeS d , ( P )€R : 



>2il 



by quadratic Hamiltonians is well known pf|. By S^, we mean the d— dimensional 
Siegel upper half space, i.e., the set of symmetric d X d matrices with positive, non- 
degenerate imaginary part [5]. The quadratic form Z describes the shape of the 
wave packet and it should be underlined that it is independent of H. In its Wigner 

representation (see eq.(J3])) ifj^'^ is a Gaufiian centered around X = y ^ ^ . This 

Wigner function is given by 



W[^](Y) = (i) e^(-\{Y-X) T G{Y-X\ 



where 



(6) G 



-^(z)^(zy 1 Q(z) + ^R(z)^(zy 1 ^(z) 



is independent of h. 

The unitary evolution, U2 (i), generated by H2 (t), acts on a squeezed state by 
translation and metaplectic action, i.e., 

(7) U 2 {t) = erMt)f(X t )M(S t ) 



1 We will, with some lack of rigor, call these states Gauftian, coherent or squeezed without 
distinction. 
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where T is the translation operator on R 2d , i.e., VT G R 2<i , —: ( ^ ^ and u G 
L 2 (R rf ) we have 

f{Y)W[u]{Z) := e^ T *- qT v)w[u]{Z) 
= W[u](Z-Y), 

where we denote by p the momentum operator and by x the position operator. The 
metaplectic operator M (F) is the quantization of a linear symplectomorphism 
on R 2d given by the symplectic matrix F. These operators form a double-valued 
unitary representations of the linear symplectomorphism of H 2d . St denotes the 
flow differential. The classical flow & H is a symplectomorphism which ensures that 
the flow differential is a symplectic matrix. Furthermore, St satisfies 

(8) ^ = JH"(X t )S u 

So — l2dx2d- 

In the prefactor of eq. (7]) we have used 

9 (t) := W (t) + fin, 

where 



W(t) := / (pTq r -H(p T , q T )) dr, 
Jo 



is the action of the classical trajectory and \x is the Maslov index of the classical 
trajectory. We have here expressed the solution of Hamilton's equations in terms 
of the canonical variables 



Pi 
(It 



X, 



Acting with J7]) on ([5]), one obtains at time t a new Gaufiian state e^ ^) ip^ up 
to a phase. The Gaufiian is centered around X t in phase space and has a quadratic 
form Zt G S^, is given explicitly by the group action [8] 

Z t = St[Z ) 

= (A t Z Q + B t ) (C t Zo + DtY 1 ,Z GS d , 

i.e., linear fractional transformation on S^. The matrix St is here written by means 
of the d x d blocks A t , B t , Ct and D t , i.e., 

V S ' " (« D, 

The quadratic form of the Wigner transform of ipg* , is explicitly given by 

(10) G t = (St l ) T G St X . 

The evolution of the approximate wave packet can be perceived as being gener- 
ated by rotation and scaling of the Gaufiian profile in R 2d . 

The difficulty of this scheme resides in the control of errors made by the approx- 
imation of U (t) in terms of U% (t). We do not dwell on the details, but just state 
the result (for proof and a thorough description see p]). The approximation can be 



2 A thorough description of the action of the evolution generated by quadratic operators on 
Gaufiians and the metaplectic representation can be found in [§]. 
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checked to all orders in ft. For this purpose, one defines the following approximating 
state 



*W(t,a:) := 



0<j<N 



where TTj ^X t , -^=, tj are polynomials in x j \fh and X t of degree smaller than or 
equal to 3j with time dependent coefficients. 

Theorem 1. (Combescure, Robert) Under the above mentioned assumptions and 
for an initial Gaufiian state tp z ° centered in the phase space representation at Xo 6 
H 2d , for every N € IN there exists C < oo such that Vft <G (0, fto] , Hq > 0, 



(11) 



u (t) i) Xa 



Z 



(t) 



N + l n . 

< C N h—te 3lt 



where < 7 < 00 is £fte Lyapunov exponent of the classical motion. 

We recall that a Lyapunov exponent is a measure of the exponential stability 
of the solutions of a differential equation upon change of initial conditions. In the 
case of classical motion, this is given by the Lyapunov exponent defined as 

(12) - ._ \r. /-I—')) 
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max 

k 



lim sup 

t — >oo 



t 



where Sk (t) are the singular values of St . The Lyapunov exponent 7 hence satisfies 
(13) IN| HS < coe^l'l 

where cq < 00 is a positive constant. We denote by ||M|| HS = ^tr (M^M) the 
Hilbert-Schmidt norm of the matrix M. Hermitian conjugation is denoted by M\ 
and transposition by M T '. 

1.3. Ehrenfest time and spreading of wave packets. It is customary, in this 
context, to define what is known as the Ehrenfest time. The latter is a time scale up 
to which the above approximation is valid. We define the Ehrenfest time, denoted 
Te (h), as the maximal time up to which 



the error 
the exact state remains localized. 



remains small, 



The latter ensures that the approximation retains a physical meaning, i.e., the above 
classical approximation makes no sense if the state does not remain localized. 
The total width of the semiclassically evolved Gaufiian is 

(14) a t (X ; h) := Ax 2 (ft, t) + A P 2 (ft, t) = Jtr(G t ), 



where we have used the variance 



Ax 2 



L 2 (H d 



2 

L 2 (R d 



and similarly for the momentum operator. 
From eq. (fT0|) one obtains that 



o* < o-q \\S t \\ ns . 
By Theorem [H there will exist a constant c\ such that 

f/(i)Vfo -e* e(t) ^ < ci\/ftie 37t . 

The time scale for which the errors are small is thus algebraic in ft if 7 = 0. In the 
generic case, the errors remain small for logarithmic times in ft. 
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With this result, the errors 
(15A(t) := (U(t)<po,s(u(t)<po)) -(u»(t)w,s(U2(t)<po)) 



L 2 (R d ) \ V / / L 2 (R d ) 

for propagating observables s can be approximated explicitly. Of particular interest 

to us is the width operator s := OpJJ' |F| 2 . One can characterize the times for 

which the error A (t) is small, e.g., A (i) = (ft Q ) for some a > 0. Again, using 
the Lyapunov inequality (fT3|) one finds that the width of the approximate state is 
bounded by 

o-t(ft) < cse 271 * 1 , c 3 >0, 

and the error remains O (h a ) up to times of order ^ ■ We may thus generically 
state that the Ehrenfest time is 

(16) r B (S oc — . 

07 

In the integrable case, implying 7 = 0, the width grows like the square of the 
Hilbert-Schmidt norm of the flow differential, i.e., at most polynomially in time. 
The error remains small for times up to ft" . The Ehrenfest time thus is algebraic 
in H. 

Our aim is to characterize the spreading of the approximate state for a more 
specific class of classical motions. 



2. Results 

In addition to the conditions imposed on symbols H above (see section II. ip , we 
assume that the Hessian H" (X t ) is T— periodic, i.e., 

(17) H"{X t ) = H"{X t+T ), VteR. 

We will also utilize the following definition. 

Definition 1 . A classical revival at a time t > is the event that the approximate 
Gaufiian given above is the initial one up to a phase factor, i.e., 

2.1. Floquet theory. According to the Floquet theorem [9], any linear differential 
equation with continuous T— periodic coefficients has a periodic solution of the 
second type, i.e., a solutions which satisfy 

f(t + T) = vf(t), ve C, Vt S R. 

In particular, the linear vector differential equation 

If - Am 

where A (i) is continuous and satisfies A (t + T) = A (t) , Vt S R, has a fundamental 
matrix of the form (see e.g. [10j) 

F t = M~ l e Lt MU (t) , M e GL (n, C) , 

where L is a diagonal matrix and U (t) is a T— periodic matrix. By definition, 
a fundamental matrix is a full rank matrix whose columns satisfy the differential 
equation, i.e., the linear combinations of the columns of the fundamental matrix 
span the full space of solutions of the differential equation. 

We will call the elements of L the Floquet exponents of the fundamental system. 

The result of this section will be summarized in the following way. 
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Proposition 1. If H" (X t ) is periodic and the Floquet exponents of St are purely 
complex then the width of a Gaufiian propagated semiclassically by U2 (t) will re- 
main unchanged at multiples of the smallest classical period T . Furthermore, if the 
classical flow is periodic then classical revivals will occur. 

Proof . Under the condition (f!7|) . Floquet theory states the existence of a funda- 
mental Floquet matrix for the linear differential eq. ((8]), 

S t = M~ 1 e At MF t , M e GL (2d, C) , 

where e At is the diagonal matrix with entries e 27rAi T ? \ i g C, F t has minimal period 
T, and we have chosen 

Fq = Sq — l2dx2d- 

One directly concludes that 

S kT = M- 1 e kTA M, VfceZ. 

for multiples k of the classical period T . 

Furthermore, the real fundamental matrix defined by S t € Sp (2d, R) has a 
unique polar decomposition, i.e., there exists [8] an orthogonal matrix Q t 6 O (2d)n 
Sp (2d, R) and a positive definite matrix Vt & Sp (2d, R) such that 

St = QtV t . 

The width of the approximate squeezed state at time t is (see eq. fl4|) ) 

a t (h) - tr((5r 1 ) T G ^r 1 ), 
and since any symplectic matrix A satisfies [8] 

A- 1 = JA T J~\ 

we can write by eq. lflO|) . 

o t (h) = ^tr (Q t VtJ T G JVjQf) = |tr (V t J T G JV?) , 

since Q t € O (2d) . 

We are hence confronted with two cases. 

• Either SkT is orthogonal, i.e., VkT = ^2dx2d', this corresponds to 
a kT (h) = hi (G ) = do (h), VkeZ. 

The orthogonality of SkT implies that it has 2d singular values 1. It is 
furthermore similar to e hTA (F k T = l2dx2d), which hence is unitary so the 
Floquet exponents are purely complex or zero. By definition, the 

Lyapunov exponent of the classical trajectories is 

7 := max lim sup , 

k [t— >oo V t ) 

where Sk (t) are the singular values of the flow differential S t , i.e., the 
eigenvalues of 

F} (M~ 1 e tA M) t M~ 1 e tA MF t = F?F t . 
Noting that F t is bounded (and periodic), we see directly that if 

max(&(Ai)) = 0, 

the classical motion is linearly stable. 
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• VkT l2dx2d, which corresponds to 

fffcr(ft) = |tr {V kT J T G JVl T ) 

> to (h) . 

since V% £ Sp (2d, R) is strictly positive definite. This corresponds to the 
case when !ft(A;) ^ 0, hence max.; (sj) > 0. The classical motion is hence 
not linearly stable in this case. 

Furthermore, if the (purely complex) Floquet exponents are rationally depen- 
dent, there exist some multiples 

Tr = UrT, n R e Z, 

of the classical period T such that the orthogonal transformation at times Tr re- 
duces to unity. Indeed, if tir is the smallest common multiple of the denominators 
of the sequence {A;} defined by the Floquet exponents, i.e., 

\iTiR e IN, Vi, 

we find 

AI~ 1 e ATR M = Af-Miagje^^jM, fijSZ, 

= ^2dx2d- 

The semiclassical approximation then is the initial Gaufiian, if it is localized at the 
initial point X . This is the case if the flow is periodic. In this case we have a 
classical revival. □ 



Remark 1. Note that these revivals are purely classical in the sense that the con- 
ditions only reflect properties of the classical motion and are hence independent of 

h. 

Remark 2. If SkT is orthogonal, the approximate Gaufiian profile will be the initial 
one at classical periods. It is just rotated and translated in phase space. 

2.2. A uniform bound. We can further characterize the approximate state in the 
case when the Floquet exponents are purely complex. 

Lemma 1. If H" (X t ) is T— periodic and the Floquet exponents of St are purely 
complex then the width a t (h) of a Gaufiian state propagated by U2 (t) satisfies the 
following uniform bound 

(18) a t (h) < Ka Q {K) , Vt G R. 

Furthermore, K := e K , where k is defined by 

k := 2T sup \\JH"(X t )\\ HS . 

te[o,T] 



Proof . We consider eq. 



^ = JH"{X t )St 



So = 1. 

Starting at some initial time kT, k G Z, St+kT satisfies the Gronwall inequality 

HWIIhs < II^Hhs e^, Vt e [0,T], 
during a classical period T where 

n = 2T sup \\JH"{X t )\\ w . 

*e[o,T] 
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Furthermore, 



<r t (h) = '^tv{S t J T G JSj) 



h 

^tr(G S?S t )< ^|tr(G )|||S«||L 



Hence 

(19) <r t (h) < ^|tr( g0 )|e K , 

if S t has purely complex Floquet exponents since in this case cr [ft) = <7kT {ft) , 
according to proposition [TJ 

We conclude with this choice that 

a t {ft) < Kctq (h) 

uniformly in time. □ 

If St is not orthogonal at the classical periods T, the latter bound takes the form 

<Tt{ft) < K<r {H)e 2max ^ h »W 

where U := are the corresponding Floquet exponents of St which is nothing 
else than a Lyapunov inequality where max, (3? (Zj)) plays the role of the Lyapunov 
exponent. Using the fact that what we are approximating is the Hilbert-Schmidt 
norm of the flow differential one can furthermore state the existence of a constant 
C4 > such that 



g 2 max;(!R(i;))|t| 



< at (h)< C4e 2max lWl ))|t|_ 



C4 

which determines the asymptotic behavior of the width in that case. 

Remark 3. We note that the equality in eq. (jTHJ) is reached if the classical period 
is null, i.e., if the Hamiltonian has a constant Hessian H" (X t ). This is satisfied by 
the harmonic oscillator. The flow differential St is then an orthogonal matrix for 
all times. It is well known that the propagation is dispersion-less in this case. 



3. Discussion 

Theorem Q] allows us to trace back the properties of the approximation to the 
exact state U (t) ip z ° . According to the discussion in [l] (see also section ll.3|l , 
this can be done with an error of order O {h a ) , a > up to times if the 
Floquet exponents have zero real part (stable classical dynamics) and up to times 
, v = maxi 3? (Zj) otherwise (unstable classical dynamics) as H \ and t f oo. 

In the stable case, the approximate state remains localized ad infinitum if this 
state remains localized between classical periods since the width of such a state is 
the same at classical periods. We wish to stress that we have defined the classical 
period as the minimal period of the Hessian H" (X t ) which is not necessarily the 
period of the classical flow & H . An example is classical motion in a one dimensional 
periodic potential Vr with energies E such that 

E > sup{V r {x)), 

In such a case H" {X t ) is periodic although the flow $^ (Ao) is not. Furthermore, 
the approximation remains localized between classical periods if 

k = 2T sup \\JH"{X t )\\ HS , 
te[o,T] 
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the exponent in eq. (fl9| . is small enough. We can hence state in that case, up to 
small errors O (h a ), localization of the exact propagation up to the Ehrenfest time 
Te cth~2. Note that the situation includes unbounded motion. 

If the flow differential St is periodic, it is clear that the shape will be the initial 
one at classical periods up to a small error. The resonance condition allows to state 
that this will also occur at some time if the flow differential is merely periodic of 
the second type and if the Floquet exponents are rationally dependent and purely 
imaginary. We may state this since the exact state remains localized and that 
our asymptotic considerations are valid as ft \ and t /* oo. The property 
extends also to the general stable case, i.e., in the case of rationally independent 
Floquet exponents. Recall that the quadratic form Z t is given by linear fractional 
transformation. In particular, 

ZkT = {AkrZo + BkT) (CkrZo + Dkr) 
where we have used the notation of section 11.21 Since SkT ~ e AfcT which is unitary 
with eigenvalues e lli ~ = e Ai ^" where ii £ R, there exist k— independent matrices 
<ii £ <C dxd such that ZkT can be viewed as the image of the vector i £ T d with 
entries under the continuous map 

e : T d ^S d 

d 

a^e 

i=0 

Furthermore, the endomorphism 

t : T d — > T d 
Lh L + i, 

where (Li + ii) is defined modulo 1, is known to be ergodic with respect to Lebesgue 
measure on the d— torus since the frequencies ii are rationally independent. In this 
notation we have Zq = e (0) = e (L (0)). From the ergodicity of r on T d we may 
state that for every e > there exist K £ Z such that \\t k (L (0)) — L (0)|| < e. It 
follows that for every e > there exist some n £ Z such that |je o t" (0) — e (0)|| < e 
implying 

Ve>0, BneZ, \\Z nT - Z a \\ < e, 

i.e., ZkT is quasi-periodic. One concludes that the initial shape of the approximated 
state will up to some small error reoccur at some multiple of the classical period if 
the Floquet exponents of the flow differential are purely imaginary. We have again 
no a priori reasons to exclude unbounded motion. 

In the case of a periodic flow, we have localization at the initial point and the 
initial profile (up to some small error) at classical periods. For bounded classical 
motion, those recurrences define revivals at periods of the classical flow that have 
already been described [TT] in the past. 

We summarize our findings in the following theorem. 

Theorem 2. Under the following assumptions 

(1) H (X) satisfies the conditions of section \l.l\ 

(2) H" (X t ) is T— periodic and v is the maximal real part of the Floquet expo- 
nents of the flow differential, solution of eq. 

and with k := 2Tsup tg [ , t] WSH" (^OIIhs where K :— e K , 

we can make the following statements up to O (h a ) , a > as h \ 0. 

// K is small enough then the approximation described by theorem Q] will hold 
up to times and the approximate width of the state will behave like e 2v \ l \ . 
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In particular, if the Floquet exponents are purely complex or zero then the 
semiclassical propagation described in theorem\]\ will hold to times %~*. The width 
of the approximated state will be bounded and recurrences will, up to small errors, 
take place at a multiple of T . 
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